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Abstract. The process of data taking, reduction and cal- 
ibration of near-infrared imaging polarimetry data taken 
with the ESO Adaptive Optics System ADONIS is de- 
scribed. The ADONIS polarimetric facility is provided by 
a rotating wire grid polarizer. Images were taken at incre- 
ments of 22.5° of polarizer rotation from to 180°, over- 
sampling the polarization curve but allowing the effects of 
photometric variations to be assessed. Several strategies 
to remove the detector signature are described. The in- 
strumental polarization was determined, by observations 
of stars of negligible polarization, to be 1.7% at J, H and 
K bands. The lack of availability of unpolarized standard 
stars in the IR, in particular which are not too bright as to 
saturate current IR detectors, is highlighted. The process 
of making polarization maps is described. Experiments at 
restoring polarimetry data, in order to reach diffraction 
limited polarization, are outlined, with particular refer- 
ence to data on the Homunculus reflection nebula around 
r] Carinae. 
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1. Introduction 

High angular resolution techniques for astronomical imag- 
ing have matured rapidly in recent years (see e.g. Beich- 
man & Ridgway 1991) and have been applied to a variety 
of Galactic and extra-galactic sources. For ground-based 
near-IR imaging, adaptive optics techniques have been 
very succesful in approaching the theo retical diffraction 
limit of the telescope (Rigaut et al. 1998 ). The systems al- 
low on-line correction for atmospheric perturbations using 
either a natural guide star (part of the object under study, 
or an unrelated star in the near v icinity) or laser guide star 
(see e.g. Lloyd-Hart et al 



observations in the near-infrared with the, in principle, 
simple provision of a polarizer in the beam. The combina- 
tion of both techniques allows information on the detailed 
polarization of an extended source, or determination of the 
individual polarization of close multiple sources. It may be 
applied to reflection nebulae for determining the position 
of embedded illuminating sources, for study of the line 
of sight geometry of dust scattering regions and for the 
orientation of magnetic fields in star forming regions or 
quasar jets. 

The extinction of interstellar dust peaks in the UV and 



declines to longer wavelengths (e.g. Mathis 199C), but the 



continuum emission from grains at typical temperatures of 
a few hundred K in regions heated by starlight increases 
strongly above 2/Ltm. In addition molecular emission and 
absorption bands are stronger above 3fim. The 1-2/im re- 
gion therefore provides an ideal window for the study of 
the close environment of dust embedded sources, such as 
regions around proto-stars or emerging young stars. For 
typical interstellar grains the low extinction in the near- 
IR enables information on the scattering properties of the 
grains, or the study of scattering regions, which have high 
optical extinction. Near-IR polarization is thus entirely 
analogous to optical polarization study but can be ex- 
tended to more embedded environments. At longer wave- 
lengths the grain emission dominates and any polariza- 
tion of the radiation is controlled by anisotropic emission 
mechanisms such as aligned non-spherical grains (Davis & 
Greenstein 1951). Examples of IR polarimetry include: de- 
tection of extended dust disks in young stellar objects (see 
e.g. Piirola et al.|l992, for observational results and Berger 



& Menard 1997, for theoretical work); dust structures in 
AGB envelopes (e.g Sahai et al. |1997| for CRL 2688); de- 
tection of dust in interstellar jets (e.g. Hodapp 1984); mag- 
netic field structure in star forming regions (e.g. Whittet et 



199^) A natural extension of ^^"^ polarization in galaxies (e.g. Jones |1997| ) and 



this technique is to achieve two-dimensional polarimetric quasars (Sitko & Yudong |199|) . Extending polarimetry to 
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the IR also brings the potential of high spatial resolution, 
both through the dependence of diffraction on wavelength 
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and the decrease in atmospheric seeing size with wave- 
length. 

In the near-IR, the dominant contribution to polariza- 
tion is therefore from scattering of radiation by grains, and 
their finite relative size requires that Mie theory must be 
used to predict the scattering properties. However the op- 
tical properties of typical interstellar grains are fairly well 
represented by models based on laboratory and observa- 
tional data (Draine & Lee 1984), so that the scattering 



properties of interstellar grains in the near-IR can be pre- 
dicted. Whilst polarization data naturally provides geo- 
metric information on the location of illuminating sources, 
the scattering efficiency with scattering angle is required 
to derive geometric information about the line of sight lo- 
cation of the scatterers (White et al. 1980 ). For high dust 
column optical depths, multiple scattering may occur and 
has then to be modelled usin g Mon te Carlo methods (c.f. 
e.g. W itt [l977t Warren-Smith |1983|; Wh itney & Hartmann 
1992|; Fischer, Henning & Yorke |1994| Code & Whitney 



1995). 



As part of a programme to study the nature of the 
dust in the Homunculus nebula around the massive star 
(or stars) rj Carinae and determine information about the 
3-D structure of the reflection nebula, near-IR imaging 
polarimetry data were obtained with the ESO ADONIS 
system. 77 Car and the Homunculus is an ideal source for 
adaptive optics since the central point source is very bright 
and the nebula is not so extended that off-axis anisoplanic- 
ity becomes an important effect. The present paper is de- 
voted to the details and subtleties of the data collection 
and removal of the instrumental signature vital to the 
derivation of a polarization map. A following paper will 
present the scientific results on the high resolution near- 
IR polarization of i] Car and the Homunculus. Sect. 2 is 
devoted to a brief description of the ADONIS instrument; 
Sect. 3 then considers the observational strategy. The fun- 
damentals of the data reduction are described in Sect. 4 
and the polarimetric calibration of the instrument in Sect. 
5. Sect. 6 exposes the different deconvolution techniques 
applied to the data and their resulting effect on the po- 
larization maps. 

2. The ADONIS adaptive optics system 

2.1. The adaptive optics system 

ADONIS is the ADaptive Optics Near INfrared System 



(see e.g. Beuzit & Hubin 1993 or Beuzit et al. 1997) sup- 
ported by ESO for common users since December 1994 
at the F/8.1 Cassegrain focus of the La Silla 3.6 m tele- 
scope. Fig |l| shows the optical layout of the ADONIS 
adaptive optics (AO) system. A tip-tilt and a 64 element 
deformable mirror corrects the distortions of the image in 
real time and a Shack-Hartmann wavefront sensor (WFS) 
provides the difference signal for the deformable mirror 
using a bright reference star close to the object. The de- 
tector for the Shack Hartmann sensor can be chosen as 



either an intensified Reticon for bright sources (m^, < 8 
mag.) or an electron bombarded CCD for fainter sources 
(8 < m„ < 13 mag., 25 to 200Hz sampling). Both detec- 
tors are sensitive in the visible wavelength region. An off- 
axis tiltable mirror allows the sky background, in a field 
of radius <30", to be chopped with the on-source image. 
The output F/45 focus delivers the image to a near-IR 
detector - either a Rockwell 256 ^ HgC dTe array (SHARP 
II for 1-2.5/Ltm, Hofmann et al. [1995[ ) or a LIR HgCdTe 
12 8^ an ti-blooming CCD (COMIC for l-b^m, Marco et 
al. 119961) . 



2.2. The camera 

The SHARP II camera was selected for the near-IR po- 
larimetric observations. This camera has a fast shutter at 
the internal cold Lyot stop, allowing integration times as 
short as 20msec. The present observations were made with 
the standard J, H, K filter set and a narrow band 2.15/im 
continuum filter, with a width of 0.017/im, and denoted 
hereafter Kc. 

2.3. The polarizer 

The polarizer, from Graseby Inc., is a wire grid of 0.25/im 
period, on a CaF2 subtrate. It is especially designed to 
work in the spectral range 1 to 9/im and has a transmis- 
sion of 83% perpendicular to the wire grid at 1.5/im. It is 
remotely rotated by the ADOCAM control system to any 
desired absolute position angle within tolerances of 0.1°. 
This polarizer, as a pre-focal instrument, is inserted into 
the beam in front of the camera; and is not cooled. How- 
ever since the polarizer is not oriented perfectly perpen- 
dicular to the optical axis there is a small image motion 
on the detector when rotating the polarizer (see Sect. 4.4 
and Fig. 5). 

3. Observational technique 

The magnification giving a pixel scale of 0.05" has been 
selected to ensure an adequate sampling of the PSF, at 
H band. The field was thus 12.8x12.8"; for the study of 
extended sources larger than the field size it is obviously 
necessary to employ several pointings and mosaic the re- 
sulting images after basic data reduction. ADONIS has a 
limit of 30" for the radial extent of the offset sky so values 
less than [30 - half detector size] (") must be employed in 
order to have unvignetted sky background frames. Special 
care has been taken in the selection of the offset sky po- 
sition to avoid any overlapping with the extended object 
observed. For all sources, object and chopped sky images 
were obtained at each position of the polarizer. A data 
cube of 256x256 spatial pixels x M frames, where M is 
the number of object and sky frames, was acquired. Table 
|l| lists the details of the ADONIS polarimetry observations 
of the science and calibration sources. 
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Sets of chopped images were obtained at nine differ- 
ent positions of the polarizer, each 22.5° apart, from 
to 180°. The minimum number of frames required to de- 
termine the linear polarization and its position angle is 3 
(spanning more than 90° in position angle). By effectively 
oversampling the polarization curve (viz. the variation of 
detected signal with polarizer rotation angle) one can at 
least hope to average out shorter term variations in at- 
mospheric transmission in order to improve the quality of 
the polarization measurement. Expressed in terms of the 
Stoke s parameters I, Q and U (see e.g. Azzam & Bashara 



1987] ), I depends on the total signal whilst Q and U de- 



pend on the difference in signals between images taken 
at polarizer angles of 0, 90, 45 and 135°. Then the linear 
polarization p is given by, p(%) = 100 x ^/{q^ + u^) 
where q=Q/I and u=U/I. The position angle of linear po- 
larization is, e{°) = 28.648 X tan-^{U/Q) (Serkowski 



1962). Determining the polarization from > twice as many 



images as necessary leads to improvement in polarization 
accuracy provided that any photometric variations are on 
timescales different from the exposure time of individual 
images at each polarizer angle. The worst case scenario 
is when photometric variations occur on a timescale simi- 
lar to the exposure times, so that the measured difference 
signals vary wildy - the polarization determined by fitting 
a cosine curve then approaches zero. The chosen expo- 
sure times per polarizer angle were in the range 1 to 50 s 
depending on the source brightness (see Table |l|). Observ- 
ing a polarized source with the polarizer at and 180° 
polarizer positions should give the same detected counts 
and is therefore a direct way to monitor the photomet- 
ric variations during the observational sequence. Column 
8 of Table |l| lists half the difference (in percentage) be- 
tween integrated counts in the star profile for the and 
180° images (i.e. rms on the mean of the and 180° signal 
values). For R Monocerotis, the semi-stellar peak of the re- 
flection nebula NGC 2261, the aperture covers the central 
extended source (full extent 8"), whilst for OH 0739-14, a 
reflection nebula around an embedded young star, an area 
10" in size was used for the statistics. 



4. Data reduction 

The data reduction applied to AO polarimetry data con- 
sists of the removal of the detector signature and sky sub- 
traction, which is common to IR imaging in general, fol- 
lowed by registration and derivation of the polarization 
parameters. 

Removal of the detector signature 

The basic data reduction steps were performed with the 
'eclipse' package (Devillard |1997 ). Flat-fields were ac- 
quired on the twilight sky at the beginning of each night 
of observation, in an exactly similar way as for the targets, 
at nine angles of the polarizer. The integration times were 



7, 10 and 20s for the J, H and K bands respectively; no 
flat field was taken with the Kc filter. The flat field images 
must first be processed to flag bad pixels, caused by either 
permanently dead pixels or ones whose sensitivity under- 
goes large fluctuation during the exposure. Two methods 
have been employed depending on the number of frames 
available in a cube: sky variation or median threshold. 

The 'sky variation' method works on a data cube, 
with preferably many planes (^20) in order to obtain re- 
liable statistics on the variations. The standard deviation 
(cr) with frame number is computed for each pixel in the 
frame. A histogram plot of the standard deviations has a 
Gaussian shape representing the response to the, assumed 
constant, sky signal. All pixels whose response is too low 
(dead) or too high (noisy), compared to a central ±ct/2 
interval, are rejected. The 'median threshold' method can 
be applied to a small number of input frames (such as 
flat field data) and detects the presence of spikes above or 
below the local mean in each individual image indepen- 
dently. If the signal is assumed to be smooth enough, bad 
pixels are found by computing the difference between the 
image and its median filtered version, and thresholding it. 
This latter method is not as stringent as using the tempo- 
ral variation, but is the only possibility when there are an 
insufficient number of images to calculate reliable statis- 
tics. Some bad pixels may however remain in the images 
after applying the bad pixel correction by either method; 
however the number is small and they can be manually 
added to the bad pixel map. Slightly different bad pixel 
maps were found for the different positions of the polar- 
izer; which could be explained by a polarization sensitivity 
of the pixels (~1%), since the NICMOS detector sensitiv- 
ity is slightly polarization dependent, or simply by the 
random variation of hot pixels. 

Once corrected for the bad pixels, the twilight flats 
were normalized, then multiple exposures were averaged 
for the same position of the polarizer to derive the flat 
field maps. The target data cubes were corrected with the 
bad pixel map derived using the 'sky variation' method 
from the background sky frames and divided by the flat- 
field to give flat-fielded, cleaned images, where the sky 
contribution is still to be subtracted. All these operations 
were performed independently for the nine positions of the 
polarizer. 

4-2. Sky subtraction 

The sky background can be bright in the IR and may also 
be polarized so it is criticical in the case of polarimetry 
to ensure that the uncertainties introduced by sky sub- 
traction are minimized. Several tests were performed to 
determine the impact of the method of sky subtraction, 
in conjunction with the bad pixel correction, on the data. 
The first method considers one sky and a bad pixel map 
for each position of the polarizer; the second method a sin- 
gle averaged sky (all polarizer positions confounded) but 
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Source 


Type 


Date 


Band 


No. 


Texp 


No. Poln. 


0-180° 










Frms. 


(ms) 


sequence 


semi-difFerence (%) 


HD 93737 


Low poln. 


1996 Mar 02 


K 


20 


50 


3 


0.76,0.65,0.22 




standard. 


1996 Mar 03 


H 


20 


50 


2 


2.32,1.01 






1996 Mar 04 


J 


20 


40 


1 


0.15 


HD 64299 


Low poln. 


1996 Mar 03 


J 


10 


5000 


1 


0.11 




star 




H 


10 


3000 


1 


1.01 








K 


10 


3000 


1 


0.89 






1996 Mar 04 


J 


3 


10000 


1 


0.05 








H 


3 


6000 


1 


0.05 








K 


5 


10000 


1 


0.65 


HD 94510 


Low poln. 


1996 Mar 04 


K, 


30 


40000 


1 


0.32 




star 














OH 0739-41 


Extended IR 


1996 Mar 02 


J 


4 


30000 


1 


0.45 




poln. source 




H 


4 


5000 


1 


0.48 








K 


4 


5000 


1 


0.20 


R Monocerotis 


Extended IR 


1996 Mar 03 


J 


10 


1000 


1 


1.26 




poln. source 




H 


20 


400 


1 


0.24 








K 


30 


100 


1 


0.18 


rj Carinae 


Polarized 


1996 Mar 02 


K 


200 


50 


4 






source 


1996 Mar 03 


H 


200 


50 


2 










H 


100 


50 


1 








1996 Mar 04 


J 


200 


50 


2 










K, 


100 


50 


2 





Table 1. List of polarization sources observed. The exposure time {T^xp) is given per frame. 



individual bad pixel maps for each position; whilst the 
third method uses the same averaged sky and bad pixel 
map for all polarizer angles. 



All three methods were tested (Ageorges 1997) and the 
results demonstrated that the largest modification of pixel 
values, and therefore photometry, comes from the bad 
pixel map used. The third method produced the largest 
discrepancies from the expected cos{26) curve, where 6 is 
the polarizer position angle. The first method is clearly 
to be preferred since the effect of any polarization of the 
sky signal on the target data is correctly removed and any 
short term variation in sky background is subtracted. 

It was found, from sky background level in the polar- 
ization calibrator data, that the sky subtraction has been 
successful to better than 1% (rms noise of 3.5 ADUs). For 
the and 180° data, a further test of the quality of the sky 
subtraction was performed: the skies have been exchanged, 
i.e. 'sky 0' has been used for the data taken at PA 180° 
and conversely. This resulted in 'photometric' variations 
less than 0.05%, thus giving us further confidence in our 
sky subtraction method. 

4-3. Photometric quality 

The photometric quality of the data can be checked in two 
different ways: either by comparing the photometry of an 
object when acquired at 0° and at 180° or by plotting 
the measured signal against the polarizer angle where a 
cos{26) form should be obtained for polarized data. The 
latter is illustrated in Fig. 0, for J band data of the NE 



lobe of the Homunculus nebula around 77 Carinae. The 
signal is plotted with time as the polarizer was rotated 
from to 180°; every ensemble of 200 points (within the 
dashed vertical lines) corresponds to frames acquired at 
the same position of the polarizer. The spread of points at 
a given polarizer angle gives a measure of the photometric 
variation. 

The images, used to create this plot, have been overex- 
posed on purpose in order to get as much signal as possible 
on the faint nebula. The central region of the images has 
thus been obtained outside the linear regime of the CCD. 
The intensity variation over this image has thus been re- 
calculated avoiding a 30x30 pixels area centered on Car. 
This is represented Fig. || together with a plot of the inten- 
sity variation over a 50 x 50 pixels area centered on a lobe 
of the nebula, away from 77 Car and thus obtained in the 
linear regime of the CCD. observed above is reduced by a 
factor of 2. Fig. 3, representing the photometric variation 
of frames acquired at and 180°, clearly illustrates the 
fact that the night of these observations was not photo- 
metric: there is a 0.3mag. extinction of the data acquired 
at 0° compared to that at 180°. 

In Fig. ^ it is clear that there is a discrepant point, 
at 157.5°, since this does not fit into the smooth cos(29) 
progression of the curve. This problem, found for every 
source observed, was attributed to a technical problem of 
unknown origin; it appears from the figure that the polar- 
izer may actually have been at an angle of 45°. All maps 
taken at this polarizer angle were ignored in the subse- 
quent derivation of polarization parameters, thus reduc- 
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ing the number of independent polarizer angles to 7 (0 
and 180° being equivalent). 

4-.4- Derivation of polarization maps 

The polarization degree for each pixel, binned pixel area 
or within an aperture was determined by fitting a cos (20) 
curve to the variation of signal with polarizer rotation 
angle 6 for the eight signal values (excluding the value 
at 157.5°). A least-squares procedure was used with lin- 
earization of the fitting function and weighting by the 
inverse square of the errors (Bcvington 1969 ). The er- 
ror on the polarization was determined from the inverted 
curvature matrix and the error on the position angle by 
the classical expression (Serkowski 1962): ag{deg.) — 
28.648(crp/p) when <Jp/p was >8 or from the error dis- 
tribution of ag /O given by Naghizadeh-Khouei & Clarke 
( |1993D when Up/p < 8. The errors on the individual points 
in the images at each polarizer rotation angle take into ac- 
count the number of images averaged, the read-out noise 
and the sky background contribution. Since the detector 
offset is not fixed per image it was necessary to boot- 
strap for the value of the sky level. A series of polariza- 
tion maps were made with increasing sky contribution at 
a fixed polarization error per pixel. The sky signal was 
adopted when it produced polarization vectors which be- 
gan to deviate from the expected centrosymmetric pattern 
(e.g. to the NE of R Mon - see Fig. ^) in the regions of 
lowest signal. Thus the polarization errors are not abso- 
lute errors. Applying a polarization error cut-off to the 
maps produces maps consistent with the expected struc- 
ture (which can also be partially checked by binning the 
data). Fig. 4 shows a typical fit to the cos{26) curve for 
a 8x8 pixels binned region of the R Monocerotis H band 
image (see Table and Fig. |). The error bars on the 
individual points arise from the photon statistics on the 
object and sky frames, with read-out noise considered. 



It was noted in Sect. 2.3 that the rotation of the po- 



larizer induces an image shift on the detector. Fig. 5 is 
an illustration of the displacement observed, for images of 
rj Carinae in Kc, while rotating the polarizer from 0° to 
180° in steps of 22.5° (see Sect. | for details on the ob- 
servation procedure). Since the PSF is variable in time, 
reproducability is not guaranteed. However the displace- 
ments were found to agree with those in Fig. 5 for different 
targets (mostly unpolarized standard stars - see Table |l|), 
and in different filters, to better than 0.5 pixel and so 
were adopted to register the images at different polarizer 
angles. 

For a point source, where only the integrated polariza- 
tion is of interest, the exact position of the source is not 
relevant provided all the signal is included in the sum- 
ming aperture. However for extended sources, such as for 
rj Carinae and the Homunculus nebula, a polarization map 
which exploits the available spatial resolution is desired. It 
is therefore extremely important to ensure that the data 



are centered on the same position for all position angles 
observed, to avoid some smearing of the information. For 
unsaturated stellar images, the centroid of the point source 
can be used as a fiducial to shift the images to a common 
centre. In the case of saturated images it proved possible 
to obtain reliable centering by using a very large aperture 
for the centroid; this is then weighted by the outer (un- 
saturated) regions of the PSF. However if the source is 
polarized, and in particular if there is polarization struc- 
ture across the point source then centroids at particular 
angles will be dependent on the source polarization. It was 
found that if the images were shifted to match the cen- 
troids at the 8 polarizer angles for the R Mon data, then a 
map with uniform, almost zero, polarization was derived, 
in contradiction to the known (aper ture) polarimetry of 
this source (e.g. Minchin et al. 1991). In such a case the 



set of image shifts, derived from unpolarized point sources 
(Fig. 5), were applied to the data and the polarization 
maps were determined. Fig. |^ shows the resulting J, H 
and K polarization vector maps superposed on logarith- 
mic intensity plots; the raw data has been binned 4x4 
pixels, i.e. 0.2". Those shifts applied are closer to reality 
than those determined by the centroid of R Mon, but good 
to within ±0.5 pixel. This might explain the difference in 
structure between our H band map and that of Close et 
al. (1997). extract of the Close et al. map. Considerable 
structure across the central (almost point) source is evi- 
dent. The cut-off of the maps is determined by the value 
of the 1(7 polarization error (4, 4 and 6% respectively for 
J, H and K). The structures seen in the J, H & K band 
maps (Fig. |^) change with wavelength, which might be an 
optical depth effect of the inclined disk. The striking dif- 
ference between the m aps in Fig. |6|and the one reproduced 
in Ageorges & Walsh (1997) comes from the calibration of 
the data. Indeed the latter were preliminary results and 
the first polarization maps derived with ADONIS. 



5. Calibration of the polarimetric data 

In order to determine the source intrinsic polarization and 
its position angle, several corrections are necessary. The 
instrument possesses an instrumental polarization which 
must be vectorially subtracted from the measured po- 
larization. The instrumental contribution is derived from 
the observation of unpolarized standards. The interstellar 
medium between the source and the observer also pos- 
sesses an intrinsic polarization which needs to be cor- 
rected. The typical ISM polarization values are <2% 
and can be neglected when observing high polarization 
sources. If the ISM polarization is not negligible, then it 
must be determined from measurements of stars in the 



neighbourhood of the source (see e.g. Vrba et al. 1976); 



alternatively the distance dependence of the ISM polar- 
ization must be determined from measurements of many 
stars. The zero point of the polarization position angle 
is checked by observing non variable polarized standards. 
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or polarized sources with reliable measurements. The lat- 
test offer an excellent check on the polarizing efficiency of 
the instrument (i.e. response to a 100% polarized source 
should be 100%). 

5.1. Sky polarization 

In the optical during dark time the sky polarization is 
typically 3-4% (Scarrott, private communication). In the 
nights of our measurements, the sky polarization has been 
found to be consistent with zero within the error bars (typ- 
ically <0.5%). Since it is the ratio of polarized intensity 
between the source and the sky that matters most and 
sinc e th e latter have carefully been subtracted (see Sec- 
tion |4.2| ) , the sky contribution has been ignored in pro- 
cessing the data. 

5.2. Instrumental calibration 

5.2.1. Choice of the polarization calibrators 

Despite extensive polarization observations, there is a dis- 
tinct lack of any such standards in the IR. The polarized 
reflection nebulae OH 0739-41 and NGC 2261 (illuminated 
by R Monocerotis) were observed because of th eir exten- 
sive I R pol arization data (Heckert & Zeilik 
et al. |l995| for OH 0739-14, Minchin et al 



Carinae had excellent spatial resolution and was retained 



(Walsh & Ageorges, 1999) 



1983 and Shure 
1991 and Close 



et al. 1997 for R Mon), although neither can be claimed 
as true, non-varying standards. 

Since the observations are achieved using an adaptive 
optics system, the polarization standard could also be used 
as a PSF calibrator. Since the correction is optimized con- 
tinuously, the resulting PSF is variable in time. Any point 
source observed as PSF calibrator needs to be close (< 
10°) to the target and be as similar as possible in terms 
of visible magnitude and spectral type, to ensure identi- 
cal correction efficiency. Owing to the lack of polarization 
standards in the infrared, the polarization calibrators were 
chosen to be as close as possible to the source and bright 
enough to be used as reference for the wavefront sensor. 
In two cases, for HD 64299 and HD 95410, which have, 
respectively a B polarization of 0.151% (Turnshek et al. 



199C) and a V polarization of 0.004% (Tinbergen 1979) 



it was assumed that the IR polarization is negligible, al- 
though no measurements exist at these wavelengths. 

In reducing the data taken on 1996 March 02 it was 
found that the derived polarization for any source (even 
OH 0739-14) was consistent with zero polarization, and, in 
addition, did not exhibit the expected shift of image cen- 
troid with polarizer angle (Fig. 5. Either the photometric 
conditions were exceptionally poor (this is not borne out 
by large discrepancies between the and 180° signal values 
- see Table |l|) or, more probably, an instrumental problem, 
such as the polarizer not rotating to the requested angle, 
was present. The polarization information was therefore 
discarded for this night. However the K band image of rj 



5.2.2. ADONIS instrumental polarization 

For the unpolarized (actually low polarization) standards, 
the integrated counts within a circular aperture includ- 
ing all the flux from the star profile (radius typically 2") 
above the sky background was measured for each angle of 
the polarizer and a cos{29) curve fitted to the data. Table 
^ lists the results. HD 93737 has a measured V band po- 
larization of 1.07% at position angle 122.4° (Mathewson & 
Ford 1970). Given the typical shape of the interstellar ex- 
tinction curve (the 'Serkowski law', see e.g. Whittet 1993 ), 
the probable values of the interstellar polarization for this 
star, assumed to have a typical Galactic interstellar ex- 
tinction, are 0.5, 0.3 and 0.2% at J, H and K respectively. 
The position angle is usually similar between the visible 



and IR (see eg. Whittet et al. 1994). For the purposes of 
computing the instrumental polarization it was assumed 
that the polarization was zero. The first two sets of data 
on HD 93737 on 1996 Mar 02 (see Table |l|) are not in- 
cluded on account of the problem with the data on that 
first night (see Sect. 5.2). In addition the first sequence 



of H band data on HD 93737 had poor photometry (see 
Table |l|) and was not considered. There is a spread in the 
values indicating typical errors of ±0.3% in linear polar- 
ization and ±15° in position angle. Given the errors the J, 
H and K values are consistent with an instrumental polar- 
ization of 1.7%. Adopted values are listed in the last row 
of the Table ^. Given that only a single measurement was 
performed at Kc, it is probably not significant that the 
instrumental polarization in this band is higher and that 
the position angle differs from the K band measurement. 

Once the instrumental polarization (intensity and an- 
gle) is determined this correction can be applied to the 
polarization maps point-by-point. Goodrich ( 1986| , in Ap- 
pendix) describes the application of the instrumental cor- 
rection. 



5.2.3. Position angle calibration 

On producing polarization maps for the Homunculus neb- 
ula around rj Carinae, it was noticed that the polarization 
vectors did not point back to the position of 77 Carinae. 
There is no reason for such a behaviour since it is known 
to be a reflection nebula. If the illumination were by an 
extended source then the offset should not be one of sim- 
ple rotation. A novel method was used to determine the 
single offset required to aligned all the polarization vec- 
tors in a centrosymmetric pattern around the position of 
rj Carinae. A least squares problem was solved to mini- 
mize the impact parameter at the position of 77 Carinae 
produced by the perpendiculars to all the polarization vec- 
tors in the Homunculus by application of a single rotation. 
A consistent value of 18±1° was found for the J, H and K 



Ageorges & Walsh: Adaptive optics polarimetry data analysis 7 



Target 


Date 


J 




H 


K 












Linear poln. (%) & PA (° 




HD 93737 


1996 Mar 02 
1996 Mar 03 






1.59, 89 


1.71, 88 




HD 64299 


1996 Mar 03 


1.51, 111 


1.99, 86 


2.16, 136 






1996 Mar 04 


1.67, 


97 


1.48, 112 


1.71, 89 




HD 94510 


1994 Mar 04 


2.12, 104 






2.05, 140 


Mean 




1.74, 105 


1.69, 96 


1.86, 104 




Adopted 




1.7, 


105 


1.7, 90 


1.7, 90 


2.0, 140 



Table 2. Polarization of low polarization stars - instrumental polarization measurement 



images. In order to verify that this was not an artifact of 
the rj Carinae nebula and the fact that the central point 
source was saturated, the 18° correction was applied to 
the polarization maps of NGC 2261. It was found that the 
vectors in the high polarization spur to the NE were well 
aligned with the direction expected for illumination by the 
peak of R Mon. Thus the calibration of the absolute po- 
sition angle can be made without reference to a polarized 
standard. 



6. Results on restoration of polarization images 

In order to measure polarization structure in the vicin- 
ity of a bright point source, it is necessary to deconvolve 
the point source response from the data frames taken at 
each position angle of the polarizer and then to form the 
polarization maps from the deconvolved images. The aim 
here is to detect polarization structure within an offset dis- 
tance of a few times the diffraction limit from the point 
source. Several different approaches to restoration have 
been attempted in order to obtain detailed information 
on the fine structure of the Homunculus nebula close to 
the central source rj Carinae. This was motivated by the 
need to detect and measure the polarization of the three 
knots found in the 0.4" vicinity of rj Car by speckle imag- 
ing in the optical (Weigelt & Ebersberger 1986| and Falcke 



et al. 1996). The polarization data for rj Car will be used 



to exemplify these experiments; the scientific conclusions 



will be reported in Walsh & Ageorges (1999). A prelim- 
inary discussion of restoration of these images, without 
considering the polarization, has been given by Ageorges 
& Walsh (|l998|). 



6.1. Image restoration trial 

Two deconvolution techniques have been applied to 
the data: Richardson-Lucy (R-L) iterative deconvolution 
(Lucy 1974, Richardson 1972|) an d blind deconvolution 
(TDAC, Jefferies & Christou |199"3| , Christou et al. |1997D . 
The major difference between these methods is related to 
the treatment of the point spread function (PSF). With 
the Richardson-Lucy method, a PSF is required a priori 
to deconvolve the data, while for blind deconvolution, the 
PSF is determined from variations in the target object 



data. The blind deconvolution method uses an initial es- 
timate, which can be a Gaussian for example. Since the 
adaptive optics PSF changes with t ime a nd is not spatially 
invariant (see e.g. Christou et al. 1998 ) , blind deconvo- 
lution should be better suited than the Richardson-Lucy 
method, which assumes a PSF constant in time. The exact 
spatial variation of the AO PSF is not known. However in 
the present case, this is a minor problem since the source 
itself (?7 Car) has been used as wavefront sensor reference 
star. Moreover with the pixel scale chosen, all the valu- 
able information in the short band data is enclosed in the 
isoplanatic angle; the spatial variation of the PSF is thus 
negligible over the area of the 77 Carinae images, which is 
not the case for the time variation. 



A comparison of the Richardson-Lucy method and 
IDAC - 'Iterative Deconvolution Algorithm in C, i.e the 
blind deconvolution algorithm used, was made using the 
Kc data on rj Car (Table 2). The aim was to test the real- 
ity of structures revealed in the near environment of the 
central star of this reflection nebula. For the R-L r estora- 
tion, the Lucy-Hook algorithm (Hook & Lucy |l994 ), in its 
software implementation under IRAF ('plucy'), was em- 
ployed. The principle is the same as for the Richardson- 
Lucy method, except that it restores in two channels, one 
for the point source and the other for the background (con- 
sidered smooth at some spatial scale) . The estimated posi- 
tion of the point source is provided and the initial guess for 
the background is flat. Kc data taken at polarizer angles 
of and 180° were restored (called K^O and Kcl80). For 
the KcO image, blind deconvolution was also performed. It 
should be noted that although the polarizer angles are ef- 
fectively identical, the Strehl ratio is not identical between 
the two data sets (K^ & K^) and is higher for Kj (27.9% 
against 22.1%). Although this could be considered an ad- 
vantage, it has a drawback since the four bumps around 
the PSF (see Fig. ^ for the appearance of the PSF) are 
more pronounced. These bumps ('waffle pattern') corre- 
spond to a null mode of the wavefront sensor as a result 
of an inadequacy in the control loop. The problem of the 
four bumps distributed symetrically around the source is 
that although they are in the PSF they do not vary; they 
are fixed in time and position and therefore not removed 
from the image as part of the PSF. There is however a way 
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Data source Polarization (%) & PA (°) 

J H K 

This work (PA unco rrected) 10.6, 77 11.1, 74 8.1, 77 

11.1, 100 8.5, 103 5.6, 102 



Minchin et al. 1991 



Table 3. JHK Polarization of R Monocerotis in an 8" aperture 



to overcome this problem, and that is by forcing them to 
be in the PSF. 

Fig. ^presents the deconvolution results obtained with 
both methods on the two separate data sets (Table |l|) and 
Fig. ^ shows the PSF derived from blind deconvolution. 
The 'plucy' deconvolved data have been restored to con- 
vergence and then convolved with a Gaussian of 3 pixels 
FWHM. The blind deconvolved data were not restored 
to convergence but limited to 1000 iterations to be com- 
parable, in terms of number of iterations, with the Lucy 
deconvolution. The resulting image seems thus more noisy 
than the Lucy deconvolved ones. Note that neither of the 
methods used succeeded in removing the 4 bumps from 
the Kc images of the first observational sequence (KjO). 

The data acquired at the polarizer position angles of 
and 180°, deconvolved with the same algorithm ('plucy') 
both show identical structures (upper row of Fig. ^. This 
example serves to illustrate the stability of the 'plucy' 
method when applied to AO data while using a reasonable 
PSF estimate. The image from the first polarizer sequence, 
polarizer angle 0°, deconvolved using ID AC is shown as 
the lower right image in Fig. ^ and is to be compared 
with the upper left image deconvolved with 'plucy'. It is 
clear that similar structures appear in both restorations 
and that there are no significant features in one restora- 
tion which do not appear in the other. The differences 
in the images are mainly due to the fact that the blind 
deconvolution has been stopped before fully resolving the 
data and the final image is thus more noisy. Moreover the 
presence of the four bumps is enhanced in this image. The 
major difficulty in this deconvolution is that these noise 
structures are convolved with extended emission from the 
Homunculus nebula. Being in the middle of the nebula, the 
flux identified on these bumps is then a convolved prod- 
uct of the waffie pattern and the extended structure of the 
nebula. It is thus very difficult for the program to isolate 
these four 'point sources' and recover properly the true 
shape of the nebula at these positions. In order to fully 
compare the different deconvolution techniques, blind de- 
convolution has been pushed to convergence for KcO (data 
set 1 & 2). The results (Fig. p|) are to be compared with 
the right hand side of Fig. ^ The structures close to 77 
Car emphasized by the two deconvolution processes, ex- 
cluding the four bumps, confer a degree of confidence in 
the scientific results which will be presented in Walsh & 



Ageorges (199£) 



6.2. Polarimetry restoration trial 

In the case of polarimetric data, the deconvolution prob- 
lem is more severe since the photometry must be preserved 
in the restored images in order to derive a polarization 
map. The Richardson-Lucy algorithm is superior to blind 
deconvolution in that it should preserve flux. Experiments 
were performed on the Kc rj Car data set, restoring each 
of the nine polarizer images with the PSF derived from 
the unpolarized standard at the same polarizer angle. The 
results were poor even when the restored image was con- 
volved with a Gaussian of 3 pixels FWHM. They illustrate 
the effect of the variable PSF and thus the difficulty to re- 
cover polarization data at high angular resolution so close 
to the star. Huge fiuctuations in the value of the polar- 
ization were seen in the vicinity of ry Car. The differing 
PSF of the unpolarized star and of ry Car (the AO cor- 
rection was much better for the rj Car images than for 
the standard star) produced restored images with large 
differences in fiux at a given pixel in the different polar- 
ization images. At present there is no known method to 
recover the true PSF from the data and conserve the flux 
through restoration. A possible (although computer inten- 
sive) solution is to determine the PSF from blind decon- 
volution and use the result for the PSF in another algo- 
rithm known to preserve the flux. This has been performed 
here: the PSF determined by blind deconvolution has been 
used both with the Richardson-Lucy and Lucy-Hook al- 
gorithms. Since the IDAC blind deconvolution algorithm 
normalised the input image at the beginning of the iter- 
ations, the final image was rescaled back to the original 
total count to allow error estimation of the polarization 
image. Polarization maps for the three methods ('IDAC 
alone and combined with R-L and 'plucy' methods) have 
been created and compared after reconvolution with a 3 
pixel Gaussian. 

From the high resolution restored images an attempt 
has been made to derive the polarization map. Fig. |l^ 
illustrates the result obtained while using the PSF deter- 
mined by the blind deconvolution with the R-L algorithm 
(30 iterations with the accelerated version), after recon- 
volution with a Gaussian of 3 pixel FWHM. The over- 
all centro-symmetric pattern of polarization observed at 
larger scale and resolution is recognisable here as well. The 
major deviation from this pattern at Aa and Ad zero (i.e. 
east-west and north-south through the image of ry Car) is 
due to the spider of the telescope. The presence of this 
feature is hard to identify on the intensity map under- 
plotted but clearly present at this position in the original 
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(undeconvolved) data. Fig. [Til is a vectorial difference be- 
tween results obtained with Lucy deconvolution and blind 
deconvolution. Special care has been taken to avoid the 
vector difference to add when the position angles were 
separated by close to 27r. Some vectors at the border of 
the noise cut-off (e.g. at Aa « —3.0") detected in the 
Lucy map but not in the other are not represented here to 
avoid confusion with the differential vectors plotted. Ma- 
jor differences can be found at approximately 0.5" from 
the center and correspond to differences in the deconvo- 
lution due to the wings of t] Carinae. At Aa = and AS 
= ± 0.3", the important difference between the two re- 
constructed polarization maps is meaningless since these 
positions correspond to the spider of the telescope and the 
data are poorly restored here. 

7. Conclusions 

The process of data acquisition and reduction for polariza- 
tion observations taken with the ESO ADONIS adaptive 
optics system has been described. Whilst certain precau- 
tions both in the observing method and in data reduction 
are required for imaging polarimetry and adaptive optics 
seperately, several other problems are presented arising 
from the combination of the two methods. 

— Since the PSF varies in time the wavefront sensor ref- 
erence star should be as similar as possible to the tar- 
get object in terms of brightness (since the achieved 
Strehl ratio strongly depends on the reference star 
magnitude) and spectral characteristic, to ensure sim- 
ilar AO correction. Since the PSF varies across the 
field of view (depending on the anisoplanatic angle) , it 
is also preferable to select the WFS reference star as 
close as possible to the target object. In practice this is 
rarely achieved when the target itself can not be used 
as WFS reference star. However a good estimate of the 
PSF provided by the reference star allows accurate de- 
convolution of the target without the introduction of 
artifacts arising from the differing PSF's. 

— Imaging polarimetry requires good photometric condi- 
tions. By oversampling the cos(2^) polarization curve 
at more than three position angles of the polarization 
analyser, an averaging over the photometric conditions 
is achieved. However depending on the time period of 
the photometric variations the averaging can result in 
zero polarization even from a substantially polarized 
source. In principle the use of an AO system should 
not compromise the photometric quality of the obser- 
vations. 

— The two polarization calibrations that are required im- 
pose the observation of an unpolarized source, to de- 
termine the instrumental polarization, and of a tar- 
get with known polarization, to calibrate the angle of 
polarization. In the IR there is a very distinct lack 
of unpolarized and polarized standards. Stars with 
known very low optical polarization are suitable as 



IR unpolarized standards since the Serkowski inter- 
stellar polarization law shows that the polarization is 
much less in the IR than the optical. However polar- 
ized standards typically have a circumstellar origin to 
their high polarization and the value in the IR cannot 
be predicted. Many of the reflection nebulae around 
Young Stellar Objects have variable polarization and 
are therefore not ideal polarized standards. 

Several strategies have been described for flat flelding and 
sky subtraction and it was shown how deviations from the 
expected cos(26') curve can give an indication of the photo- 
metric conditions at the time of observation and allow any 
discrepant polarizer angles to be discarded as was found 
for the ADONIS polarizer at PA 157.5°. The instrumen- 
tal polarization for ADONIS was determined at 1.7% over 
the J, H and K range. Polarization maps have been succes- 
fuUy produced for the reflection nebula around ij Carinae 
(the Homunculus). By using the PSF's determined from 
blind deconvolution at the same polarizer angles as the 
data, it has been shown that polarization structure can 
be revealed as close as two times the diffraction limit to 
a point source. The interpretation of the ADONIS AO 
polarization results on rj Carinae will be presented in a 
forthcoming paper (Walsh & Ageorges 1999). 
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Fig. 1. Optical layout of the ESO ADONIS adaptive optics system. The polarimeter, as a prefocal instrument, is 
installed at the entrance window of the dewar (IR focus f/45). 
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Fig. 2. The photometric variation of the basic data is illustrated by the time sequence of measured counts in a region 
of the Homunculus taken in J band. Every 200 frames, the polarizer has been rotated by 22.5° and the vertical dashed 
lines indicate the change of polarizer position angle. The width of this curve is characteristic of the photometric 
variations. The discrepant point at 157.5° is attributable to an instrumental problem (see text). Left: global intensity 
variation over the full data frame but excluding a 30x30 pixel box centered on r] Carinae; right: variation over a 50 x 
50 pixel area of the nebula, far from the saturated center of the image. 
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Fig. 3. Photometric variation between data taken at PA and 180° is shown for the full frame (excluding the central 
source, i.e. ij Car itself) at left, and for a small area centered on the Homunculus (right). 
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Fig. 4. A typical fit of the observed signal as a function of polarizer rotator angle by pcos{29) for the summed counts 
in an aperture over the H band image of R Monocerotis (see Table |l| and Figure The derived value of linear 
polarization and position angle is shown by the bold line. The point at 157.5° was not considered in the fit. 
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Fig. 5. Illustration of the variation of the position of the centroid of an image on the detector while rotating the 
polarizer from to 180° in steps of 22.5°. The target was r] Car observed in Kc band and the shifts are clockwise with 
increasing rotator angle. 
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Fig. 6. J, H and K polarization maps of R Monocerotis and the brightest regions of the reflection nebulae NGC 2261. 
The contour maps shows the logarithm of the signal (polarized + unpolarized) . The size of the polarization vectors is 
indicated and the orientation is north to the top and east left. 
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Fig. 7. Comparison of the results of two deconvolution techniques for Kc data of rj Carinae and the Honiunculus nebula. 
On the top row are shown data from the second polarizer sequence, at (left) and 180° (right), both deconvolved with 
'plucy' and convolved with a Gaussian of 3 pixels FWHM. The bottom left image is identical to the top left one but 
the data restored is from the first polarizer sequence (0° image). The result of blind deconvolution is presented in the 
bottom right hand image for the first data set - - at 0° polarizer angle. 
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Estimated PSF IDAC determined PSF 




Fig. 8. PSF estimated (left) and obtained (right) from an IDAC blind deconvolution of the Kc image of rj Car (see 
Table 2 for observational details). The PSF estimate is the shift-and-add of the point som-ce (HD 94510) observed 
shortly before the target data with the Kc filter. It has been used as first guess for blind deconvolution. 
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Fig. 9. Kc data of ?/ Car. at 0° for the two different data sets acquired. These images have been obtained after 
convergence of the blind deconvolution algorithm and reconvolution with a Gaussian of 3 pixels FWHM. This has to 
be compared with the right hand pictures of figure top and bottom respectively, which are the equivalent results 
for 'plucy' deconvolution. 
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Fig. 10. Polarization map of Kc deconvolved data of rj Car. overplotted on the high resolution intensity map. For 
clarity, the polarization vectors have been calculated on pixels binned 4 by 4 (0.2x0.2"). The intensity map is at the 
nominal resolution of 0.1" 
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Fig. 11. Differential polarization map of Kc deconvolved data of rj Carinae. It represents the difference between 
the polarization map derived after pure blind deconvolution and R-L deconvolution using the PSF derived by blind 
deconvolution. The errors are clearly under 5% except at the border and thus low signal_to_noise level of the nebula 
and at the position of the telescope spider (Aa and AS = 0). Compared to the rest of the map some bigger differences 
can also be found at « 0.5" from rj Car and is attributed to the effect of the PSF wings in the deconvolution process. 



